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PROCEDE D'EXAMEN RADIOLOGIQTJE D'UN OB JET 
DESCRIPTION 



Cette invention concerne un precede 
d'examen radiologique multi-ehergie d'un objet. 

Lea procedes radiologiques consistent a 
fair, traverser un objet a etudier par un rayonnement 
pour en deduire la repartition de differentes 
categories de materiaux, absorbant le rayonnement de 
fasons differentes, dans cet objet. Une ' application 
tris courante est 1 'osteodensitomStrie, ou on recherche 
la masse et la densite des tissus osseux dans un 
patient en distinguant ces tissus des tissus' mous. 

11 est usuel d'utiliser un spectre de 
rayonnement large et a le diviser en bandes mesurees 
separement par des canaux de mesure respectifs. Comma 
les coefficients d'absorption ou d-attenuation du 
rayonnement par 1-une ou !■ autre des categories de 
tissu sont different* pour chacune des bandes le 
probleme theorigue se ramene a . la solution d'un systeme 
d- equations dont le nombre est egal a celui des bandes 
de mesure et qui comprennent chacune deux inconnues 
des epaisseurs ou les masses traverses de tissua mous 
et osseux). La solution du probleme devient possible 
apres une calibration obtenue en faisant traverser par 
le rayonnement divers etalons dotes d- epaisseurs 
connues de materiaux " aux proprietes d'absorption 
analogues a celles des materiaux de mesure, et 
notanunent le plexiglas et 1 -hydroxyapatite pour simuler 
les tissus mous et les tissus osseux. On peut alors 
calculer lee parametres mathematics d'un modele 
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reliant las mesures d - attenuation aux epaisseurs des 
matSriaux. 

On notera qu'en rgalitg le ^ 
contxent trois categories principals de tissus : les 
tissus osseux, les tissus maigres et les tissus gras 
maxs gu-on n-en considere que deux en general a cause 
de difficulty a distinguer ces trois categories dans 
les mesures, de sorte qu-on confond volontairement les 
txssus maigres et gras. D'autres precedes sont ensuite 
appliques pour distinguer leurs proportions dans les 
tissus mous. 

La largeur du spectre permet de disposer 
d'un nombre beaucoup plus grand de bandes de mesures 
lu-zl ne serait necessaire pour obtenir une solution 
et de les employer toutes. pour obtenir des resultats 
Plus precis en exploitant tout. 1 ■ information 
d- absorption obtenue. Dans 1- article -Measurement of 
bone mineral using a multiple -energy x - ray 
absorptiometry", par j. Svanpalmer, R. Kullenberg t 
Hansson, Phys. Med. Biol, Vol 43, 1997: PP 379-387, ou on 
dispose de 23 bandes de mesures et ou on considere les^ 
trois categories de tissus, il est propose de combiner 
les groupes de mesures trois par trois de toutes les 
facons possibles pour obtenir . 1771 (23x22x21/2x3) 
systemes de trois equations .a trois inconnues qui 
donnent autant de groupes de resultats. 11 convient 
alors d'appliquer un critere de choix. Les auteurs 
conseillent de choisir comme resultat veridique celui 
qui est a la valeur moyenne, ou ft la valeur mediane 
Pour le parametre le plus important, qui peut etre la 
masse osseuse travers€e. 
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Contrairement a ce qu'on a consent plus 
haut, les auteurs considered les trois categories de 
txssus. cela ne met pas en question la validite de leur 
methode, i condition d-exposer le patient a une 
mtensiti d- irradiation beaucoup P l ua forte pour 
ramener les incertitudes des mesures a des proportions 
acceptables. II se rait d-ailleurs possible de modifier 
sxmplement leur methode pour 1-appXiquer a des mesures 
sur deux tissue, ou plus generalement deux categories 
de materiaux, seulement. Toutefois, il 'subsiste le 
defaut que cette methode ne tient pas compte du bruit 
aur les mesures et qu'elle conduit done a des resultats 
bruit^s . 

L' invention a pour objet un 
perfectionnement de tel. precedes par combinaison de 
resultats nombreux et comprend une amelioration du 
critere de choix des resultats. Plus precisement, elle 
conceme sous sa forme la plus generale un precede 
d'examen radiologique d'un objet ou sont considers au 
moins deux categories de materiaux, comprenant - 
l'emploi d'un rayonnement a spectre large ; des mesures 
du rayonnement par bandes du spectre ; des expressions 
IM) d'epaisseurs ou de masses des deux categories de • 
materiaux traverses par le ray0 nnement, les 
25 expressions (M) £t ant des fonctions d'au moins deux 
des mesures ( me s k ) et de coefficient (A) ■ et 
^application d-un critere de choix parmi ' les 
expressions (M) pour en dgduire ^ ^^.^ 
Wfinale) jugee veridique ., caracterise en ce que le 
critere de choix comprend un combinaison (f) des 
expressions avec des facteurs de pondSration (a) , et un ■ 
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calcul de facteurs de ponderation de maniere que la 
combinaison ait un bruit minimal (variance minimale en 
langage mathimatique) calcule d'apres lee bruits sur 
les mesures (variances sur les mesures) . 

L ' inVention sera nontenant decrite en 
reference aux figures : 

- la figure 1 est une vue d'un spectre, 

la figure 2 est une repartition de 

resultats, 

" St la figure 3 est un prganigramme 
resumant 1' invention, auquel on pourra se reporter dans 
toute la description qui suit. 

L- attenuation du rayonnement peut Stre 
exprimee par une fonction des epaisseurs traverses de 
15 chacun des matSriaux d- indices x et y, ou de. leurs 
masses M (densitee par unite de surface) dans la 
direction du rayonnement, Le spectre de mesures de la 
figure 1 est divise en N bandes notees generalement par 
les indices i et j . Les attenuations varieront dans 
20 chacune des bandes en raison de coefficients 
d- absorption variables pour les deux materiaux. si nous 
appelons me Si ou mes, les mesures pour une bande 
d'energie i ou j , i es masses traversees M x et My 
pourront chacune etre exprimeea aussi, par exemple, par 

25 \* f °™ le generale 

M = Al+ A 2 . mes i+ A 3 .mes j+ A, . mes » + A S . mes ) +As . me Sl .mes, . . i es 
. mesures considers dans cet exemple etant des meusures 
d' attenuation, on aura pour chaque canal de mesure i 
(correspondant a une bande du spectre)" la' relation me Sl 

3 0 = i n { noi/ ni) oa no . estle nombre ^ phatQns arr . vant 

l'objet et ni celui des photons ayant traverse l'objet. 
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Comme la non-linearite des fonctions et My en 

fonction des mesures est faible en pratique, on pent 
e-accommoder de cette fonction polynomial du deuxieme 
degre qui comprend six coefficients Al a A 6 . L e degre du 
P olyn6me peut etre ajuste en fonction du probleme , Par 
exemple, pour i <analyse d , objets ^ 

materiaux ayant des numeros atomiques plus eleves que 
elm des tissus biologiques, comme le pontile non 
destructif pour l-examen d'un objet metallique. les 
mesures considers dans cet exemple Stent des measures 
d' attenuation, on aura pour cheque canal de mesure i 
(correspondent a une bande du spectre) la relation me Si 
- m ("<%.) 0 * no . estle nombre ^ phatoas arr . vaiit ^ 

1' objet et ni celui des photons ayant traverse 1- objet. 

Ces coefficients peuvent etre trouves dans 

une etape d-etalonnage a travers des etalons/ appeles 

parfo 1S fantomes ou cales dans 1'art, et qui consistent 
en des pieces formees d'epaisseurs connues et 
dxfferentes entre elles de materiaux simulant, par 
leurs propriety d ■ attenuation, les materiaux de 
1 'objet qui sera ef f ectivement a mesurer. Chacun des 
etalons est done sounds a une irradiation par le 
rayonnement pendant une longue duree qui permet de 
reduire 1 ■ influence du bruit sur les mesures. Le 
spectre mesure pour chacun des etalons donne encore N 
mesures resultant de la decomposition . du spectre en 
autant de bandes. En combinant maintenant ' deux series 
de mesures i et j prises pour deux bandes et pour 
chacun des etalons, on recherche les coefficients A 
pour ajuster les fonctions M x et M y aux mesures. Dans le 
cas present, ou il y . six coefficients A pour chacune 
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des deux fonctions, et' ou on considere les combinaisons 
de deux bandes de mesure, les mesures pour chacune des 
bandes devront, pour foumir une solution unique 
porter sur six etalons. Des etalons plus nombreux 
5 pourront aussi Stre utilises pour ameliorer la 
predion sur la solution. Une minimisation de fonction 
d'erreur sera alors appliguee. 

Cette determination des coefficients A est 
repetee pour diverses oombinaisons de . paires de 
10 mesures. II en itait encore ainsi dans !• article 
anterieur mentionne plus haut ; cependanti on a 
constate qu-il etait inutile d-effectuer toutes les 
combinaisons, au no.bre.de ^zl>, pour exploiter 

complement les mesures et que ( N - 1} combinaison 
IS etaient en realite suffiaantes pour recueillir toute 
1 ' information. 

Une facon preferee de proceder consiste a 
choisir au depart la bande de mesures qui est la moins 
bruitee (par exemple celle qui a le signal n - le 
10 nombre de photons a la reception- le plus important) et 
a 1'associer successivement a chacune des autres bandes 
de mesures pour les combinaisons. On obtient finalement 
(N-l) ^estimations des deux fonctions M* et My, gu-on 
note M V M 2 ,...,MN~1 pour chacune de ces deux fonctions. 
5 A ce stade du precede, les fonctions M x et 

M„ representatives des longueurs traversees de deux 
materiaux representatif s des tissus osseux et des 
tissus mous, pourront etre converges en fonction Mv 
et M„ representatives des longueurs equivalentes 
3 traversees des tissus osseux, des tissus maigres et des 
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txsaus gras en confcinant lineairement * Bt ^ de trois 
f aeons differentes determines par 1 ' experience Ce 
precede de conversion est independant de 1- invention, 
dSja connu dans Hart. 

Comme il n-y a pas de raison de preferer 
l'une ou !■ autre de ces estimations A. un critere de 
choix doit etre applique pour obtenir !■ estimation 
Mfinale qui sera jugee veridigue. Dans 1'article 
anterieur, une dee expressions obtenues etait 
direct ement select ionnee d'apres .un critere de 
• classement (l e resultat median) ou de moyenne des 
valeurs prises par les expressions pour un des 
r6sultats. Dans 1' invention, les expressions M seront 
combines; par exemple lineairement, d.'apres la formule 
& finale = (al Mi) + (a2 M 2 ) + ... + (a N-i Mjs _ 1} tout 
en minimisant le bruit ; les coefficients a 1# etc. ont 
une somme egale a 1' unite (al + a2 + _..+ a N-l = i) . 

Pour chague canal de mesure, le bruit susr 
le nombre de photons suit une loi statistigue de 
Poisson, dont le resultat est independant pour chacun 
des canaux. La matrice de covariance des N-l resultats 
peut Stre exprimee d'apres la formule 

r -T d1 * 1 d ** j 1 

. >J *=i dmes k dmes k N k ' - 

La variance sur la combinaison lineaire 
dormant U finale s'exprime par la formule f =(a 1 _,k H . 
i) .rSl^, . cette quantity f atteint ^ valeur 
optimale quand sa derives selon toutes les variables 



B 14155. 3/JCI 



AUG-27-2003 10:33 



21S 579 6073 



91* 



P. 08 



8 



est nulle, c-est-a-dire que 1' influence du bruit est 

'dffda^O 
df/da 2 =0 



minimisee, eoit : < 
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La resolution numerique de ce systeme 
foumit les coefficients al, a2, etc et M finale, 
c'est-a-dire les masses traversees des deux categories 
de matSriaux. 

Ce sont ces operations qui sont accotnplies 
dans I- invention ; il est important de remarquer que 
les mesures (mes k ) employSes dans la formule donnant r f 
et ensuite sont les mesures effectuees a travers 
1' 6b jet a etudier lui-meme, mais pas celles qui ont fits 
faites a 1'etalonnage pour determiner les coefficients 
A des f onctions M . 

L- invention permet d'echapper au dilemme de 
precedes anterieurs dans lesquels on renoncait a une 
partae des energies de mesure ou au contraire on 
acceptait des mesures peu distinctes : elle utilise en 
effet tout le spectre, mais en le divisant en bandes 
assez nombreuses pour que la mesure de chacune puisse 
Stre comparee de facon utile a d-autres mesures, faites 
sur des bandes eloignees. Elle permet ainsi, entre 
autres avantages, d'etudier aussi bien des organismes 
maigres que gras . ' 

Enfin, 1- invention pent etre generalise a 
un no^ re de materiau* plus grand que deux, ce qui peut 
etre rnteressant notamment pour les precedes d-imagerie 
a produit de contraste, ou trois variables doivent Stre 
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considerees ou pour contrSle dea bagagea ■ (recherche 

d-explosifa) , et elle pent etre appliquee en 

considerant des combinaiaons autrea que lineaires des 
f onctions M . 

Dea application de 1 ' invention sont : 
. 1 'ost€odensitomitrie 

- obtention de la densite de masse 

osseuse 

obtention de la .composition 
corporelle (masse grasse, masse maigre) 
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. le controle agro-alimentaire, par exemple 
la detection de fragments d'oa dans de la yiande ou la 
detection de morceaux de verres dans des plats 
cuisinis, 

- le controle des bagages . : recherche 
d'explosifs, de produits illicites (armes, aliments, 
drogue™) 
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REVINDICATIONS 



1) Procede d'examen radiologique d'un objet 
ou sont considerees au moins deux categories de 
materiaux, comprenant : l-emploi d'un rayonnement a 
spectre large ,- des meaures du rayonnement par bandes 
du spectre ; des expressions {M) d'epaisseurs ou de 
masses des deux categories de materiaux traverseee par 
le rayonnement, les expressions (M) •• etant des 
fonctions d'au moins deux des mesures (mes k ) et de 
coefficients (A) , et 1 - application d'un critere de 
choix parmi les expressions (M) pour en deduire une 
expression (j* finale) jugee vSridique ,- caracterise en 
ce que le critere de choix comprend une combinaison (f) 
des expressions avec des facteurs de ponderation (a), 
et un calcul des facteurs de ponderation de maniere que 
la combinaison ait un bruit minimal calcule d'apres un 
bruit sur les mesures. 

2) Procede selon la revendication 1, 
caracterise en ce que la combinaison (f) des 
expressions est lineaire. 

3) Procede selon la revendication 1 
caracterise en ce que la variation de la combinaison 
est calcul.ee avec une matrice de covariance <r f ) des 

25 expressions {M). 

4) Procede selon la revendication 1, 
caracterise en ce que les expressions [&) sont en 
nombre egal a celui des bandes, moins un, et etablies 
toujours avec une des bandes Cio) et, respect ivement , 

3 0 chacune des autres bandes. 
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5) Procede selon la revendication l, 
caracterise en ce que les coefficients ( A ) sent 
determines a une etape prealable de calibration. 

6) Procede selon la revendication 1, 
5 caracterise en ce qu'il est applique i 

1 1 osteodensitometrie . 

7) Procede selon la revendication l, 
caracterise en ce qu'il est applique a des controles 
agro-alimentaires . 

10 8) Procede selon la revendication l, 

caracterise en ce qu'il est applique au'contrSle de 
bagages . 
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L« invention concerne un procSdS d' examen 
radiologic^ d'un objet oil sent considered deux 
categories de materiaux, comprenant : l^emploi d'un 
rayonnement a spectre large ; des mesures du 
rayonnement par bandes du spectre ; des expressions 
[A) d'epaisseurs ou de masses des deux categories de 
materiaux traversers par le rayonnement, les 
expressions (M) stant des fonctions d-au moine deux 
des mesures (mee k ) et de coefficients (A) ; et 
l'application d'un critere de cnoix parmi les 
expressions (M) pour en dgduire m expreBa±on 
(Mfinale) jugee veridique ; caracterise en ce que le 
critere de choix comprend une combinaison (f.) des 
expressions avec des facteurs de ponderation (a) , et un 
calcul des facteurs de ponderation de maniire que la 
combinaison ait une variation minimale d'apres des 
variations des mesures. 

Figure 3. 
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MESURES A TRAVERS DES ETALONS 
POUR UN SPECTRE A N BANDES (ij) 



CALCUL DES COEFFICIENTS (A) 
DES FONCTIONS DONNANT LES EPAISSEURS 
TRAVERSEES (CONNUES) D'APRES 
LES MESURES (mes; .mesj ) 



MESURES A TRAVERS L'OBJET A ETUDIER 



ESTIMATIONS (M) DES EPAISSEURS 
TRAVERSEES EN APPLIQUANT 
LES FONCTIONS D'ETALONNAGE 



COMBINAISON (f) DES ESTIMATIONS 
AVEC DES FACTEURS DE PONDERATION (a) 



AJUSTEMENT DES FACTEURS DE 
PONDERATION POUR MINIMISER 
UNE VARIATION DE COMBINAISON 
EN FONCTION DES MESURES 



ESTIMATION FINALE DES EPAISSEURS 
EN APPLIQUANT LA COMBINAISON (f) 
ET LES FACTEURS AJUSTES 



FIG. 3 
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